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Abstract 

We discuss the resummation of next-to-leading logarithms (NLL) 
in heavy quark production near threshold. Results are presented for 
top quark production at the Fermilab Tevatron. 

1 Introduction 

There has been a lot of interest recently in the resummation of soft gluon 
radiation near threshold for heavy quark production in hadronic collisions 
[1-4]. This resummation is a direct consequence of the factorization of short- 
from long-distance effects in QCD scattering [5]. In the perturbative expan- 
sion for the heavy quark cross section one finds logarithmic terms that can 
be resummed to all orders of perturbation theory. These logarithms come 
from distributions that are singular for z = 1, where z = Q 2 /s, with Q 2 and 
s the invariant mass squared of the produced heavy quark pair and the par- 
tons in the incoming hadrons, respectively. One of the interesting features 
of singular distributions in QCD-induced cross sections is their sensitive de- 
pendence on the color exchange in the hard scattering. This feature as well 
as the presence of final state interactions are the main complications rela- 
tive to the Drell-Yan process [6]. We will show that these effects contribute 
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only at NLL. As we discuss below, resummation is conveniently carried out 
in terms of moments, and logarithmic dependence on the moment variable 
exponentiates. The moments have to be inverted to derive the physical cross 
section. Previous formalisms [1-4] resum only leading logs (and some NLL) 
and there has been a debate as to which subleading logs to keep. In our anal- 
ysis (in moment space) we have for the first time resummed all the NLL so 
we have effectively pushed the debate in [1-3] to the next-to-next-to-leading 
logarithmic level. In Section 4 we invert moments using the method in Ref. 
1. 

2 Resummation for QCD hard scattering 

For the hadronic process ^i(pi) + h 2 (p2) — > QQ + X the cross section a, 
can be written as a convolution of parton distribution functions 0, with the 
partonic hard scattering a, 

V hl h 2 ^QQ = E <l>f/hi ® Vff^QQ ® <i>f/h2 ■ I 2 ' 1 ) 

f=q,Q,9 

This convolution becomes a product in terms of moments, 

^m»W= E h/ hl (N)tff->QQ(N)t f/h2 (N), (2.2) 
f=q,q,g 

where a{N) = ft dr r N ~ l cr(r), a(N) = ft dz z N ~ l a(z), and 0(iV) = 
Jq dx x N_1 4>{x). Here r = Q 2 /S where S is the invariant mass squared of 
the incoming hadrons. Singular distributions in z give rise under moments 
to logarithmic iV-dependence, as in 

oc ln m+1 N. (2.3) 

+ 

Thus, the large- iV behavior is a diagnostic for singular distributions in 1 — z. 
We can derive the resummation of soft gluon contributions by reexpressing 
moments of the cross section (with hi — f, h 2 — f) in a form in which 
collinear gluons are factorized into alternate parton distributions ip, and soft 
gluons into a function Sji, 



f 

Jo 



dz z 



N-l 



ln m M 



1-z 



0-//-QQ = i/>f/f ® hjhi ® Sji ® ipf/f. 



(2.4) 
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This factorization takes into account the color exchange at the hard scat- 
terings, hi and hj in the amplitude and its complex conjugate, respectively, 
where I and J are color tensor indices. 

From Eqs. (2.1) and (2.4) we find a(N) = [ip(N)/(j)(N)} 2 h} S n (N) hi. 
The ratio of the universal ip and has been analyzed in Drell-Yan production 

!• 

The soft matrix S ji(Q/(N/i)) requires renormalization: sf) = {Z\)jbSba 
xZs,ai, where Zg,ji is a matrix of renormalization constants @. Hence Sji 
satisfies the renormalization group equation 



fi- 



dS 



.ii 



d\i 



—(T s )j B Sbi — Sja(Ts)ai 



(2.5) 



In a minimal subtraction scheme, and with e = 4 — n, the anomalous dimen- 
sion matrix is 

n f) 

(2.6) 



Ts(g) = --g-ReSe-^oZs^e) 



At the level of leading logarithms of N in Sji, and therefore at NLL of N in 
the cross section as a whole, we choose a color basis in which the anomalous 
dimension matrix is diagonal, with eigenvalues A/ for each basis color tensor 
labelled by /. Then, the solution to Eq. (2.5) is 



Sji ( jj, a s (Q 2 ] 



S 



ji 



x exp 



1, a x 



Qf 

N_ 
dp 



^[Xi(a s (^)) + A}M/i 2 ))] 

Q/N [I 



.(2.7) 



Thus, we find for the partonic cross section [8,9] 
%->qq(A0 = A'e E ^ N ^hj(l,a s (Q 2 )) sji, 



a. 



hi(l,a s (Q' 



(2.8) 

where A' is an overall constant and 9 is the center-of-mass scattering angle. 
The exponent is 



Eji(N,9,Q 2 



E DY (N,Q 2 )- / dz 



yN-l 



1 r 



1-2 



+ g\ J ^[a s {{l-zfQ 2 )A] , (2-9) 
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where Ejyy is the Drell-Yan exponent, and 



gi I) [a s ,8] = -\ I [a s ,6} + ^ C F , A . 

7T 



a. 



(2-10) 



3 The Ts matrices for qq —> QQ and gg — > QQ 

We now give explicit results for the anomalous dimension matrices Ts- We 
begin with quark- ant iquark annihilation, q(p a ) + q(pb) Q(Pi) + Q(P2)- 
The explicit calculation is given in ||. Here we give the results for the 
anomalous dimension matrix in a color tensor basis consisting of singlet and 
octet exchange in the s channel, 



(3.1) 



The anomalous dimension matrix (shifted by Cj?a s /7r as in Eq. (2.10)) in 
this color basis and in an axial gauge A = is [f| 



r,S,22 



Q> 

-C F {L p + \ + m), 

IT 

2<y s /«a C F 
— ln W' r5 < 12 = ^ 

a. 



r 



5,21, 



7T 



c h 



41n(^)-L,-l 



Til 



Ca 
2 



-In 



2 2 

u\m s 
'i 



+ Lr + ixi 



(3.2) 



where s = (p a + Pb) 2 , h = (p a — Pi) 2 — m 2 , u\ = (pb — Pi) 2 — rn 2 , and m is 
the heavy quark mass. Also 

1 - 2m 2 /s ( I- (3 



In 



1 + 13 



+ 7U 



(3.3) 



where (3 = yl — Am 2 /s. T$ is diagonalized in this singlet-octet basis for 
arbitrary (3 when the scattering angle is 9 = 90° (where U\ = £i). It is also 
diagonalized at /3 = 0. 

Next we discuss gluon fusion, g(p a ) + g(Pb) —> Q(pi) +Q(P2)- We choose 
the following color basis: 



ci = S ab 5 21 , c 2 = d abc T^, c 3 = if ahc T^. 



(3.4) 
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Then in the A = gauge the anomalous dimension matrix (shifted by 
C A a s /TT as in Eq. (2.10)) is @ 



r Sl u = -[-c F (Lp + i)-c A <Ki] 

7T 

1^5,21 = 0, F$A2 = 0, 



r.csi = 2 — In 



7T V tl 

C A 



Ts,22 = {— Cf(L/3 + 1) 

7T 



+ 2 



In ( ) Lp - nl 



m?s 



N 2 -A 1 

1^5,32 — — ^ — r 5j31 , ^^13 — -is^i, 

Ts,23 = —T-Fs,31, F S ,33 = ^S,22- (3.5) 

We note that the matrix is diagonalized at 9 = 90° and also at /3 = 0. 

We have checked that the one-loop expansion of our results for both 
channels are consistent with [10]. 



4 Results for top quark production at the 
Fermilab Tevatron 

We are now interested in the magnitude of the NLL terms, particularly the 
#3 contribution, relative to the LL results of [[TJ. Therefore we use a similar 
cutoff scheme, modifying the definitions of the exponent in Eq. (2.9) to work 
directly in momentum space. This method uses the correspondence between 
logarithms of the moment variable N and logarithmic terms in the momen- 
tum space variable 1 — z. We work at 9 = 90° to avoid the diagonalization 
of the anomalous dimension matrices. 



The resummed partonic cross section is |TT 



<Jab{s, 17l 2 )=- d Si /aft hr"? \ , (4.1) 



,/scut \2m 2 J ds 4 

where ab = qq or gg and S4 = 2m 2 (1 — z) at threshold, dxrj^ jds^ is the 
differential of the Born cross section JTJ. After replacing z N ~ 1 - 1 in Eq. 
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Fig. 1. The resummed top quark production cross section at 6 — 90° in the 
MS scheme. We show the resummed results without the g 3 terms (lower solid 
line) and the resummed results with all the NLL terms (Q = to, upper solid 
line; Q = 2m, lower dashed line; Q = to/2, upper dashed line). 
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(2.9) by —1 and introducing uj = 1 — z, the exponential function for the qq 



channel is 11 



f qq (^) = exp[^ Y + E qq -(\ octct )] , (4.2) 

where E^ is the Drell-Yan part. The color- dependent g% contribution in 
Eq. (2.9) leads to 



Egg (Ai) 




+ A* 



(4.3) 

where i denotes singlet or octet. Since our calculation is not done in moment 
space, the uj integral is cut off at uj$ = s^/2m 2 . Because the running coupling 
constant diverges when u 2 Q 2 /A 2 ~ 1, the minimum cutoff in eq. (4.3) is 
-Scut — ■S4. m i I1 ~ 2m 2 A/ Q, where A is the QCD scale parameter. In general 
we choose a larger value for the cutoff consistent with the sum of the first 
few terms in the perturbative expansion; here we cut off the ^-integration at 
z = l- 10A/Q. 

The treatment of the gluon-gluon channel in the MS scheme is very similar 
but now we have three distinct color structures. However, only two of them 
are independent. We define f gg for each eigenvalue so that f g9t i = exp[2^ Y + 
Egg(Xi)}, in complete analogy with the relations for the qq channel. For 
details see Ref. 11. 

In Fig. 1 we show results for the resummed top quark cross section at 
the Fermilab Tevatron with = 1.8 TeV (with and without the NLL g 3 
terms) versus the top quark mass. We used the MRSD— ' parton densities 
[12]. We also show the variation of the cross section with the factorization 
scale. Then for m = 175 GeV/c 2 and Q = m the total ti resummed NLL 
MS cross section at 9 = 90° is 5.3 pb. Without the NLL g s contribution the 
cross section is 3.7 pb. Thus the NLL g 3 terms enhance the total MS cross 
section considerably. 



5 Conclusions 

We have given explicit results for the resummation of next-to-leading log- 
arithms in top quark production. We have shown that the NLL terms are 
numerically significant. Our methods have been applied to 6-quark produc- 



tion at HERA-B, with similar conclusions [11|, and they can also be applied 
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to jet production We have recently completed the evaluation of the 

anomalous dimension matrices for qq — > gg, qg — > qg, and the more compli- 
cated case of gg — > gg, all relevant in jet production [O . 
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Questions 

L. Trentadue, Parma: 

You have shown some numerical outputs of your calculation. You have 
used the method of Laenen, Smith, and van Neerven to cut-off the z — > 1 
limit. Could not the results be affected by this by hand limiting of the final 
phase space? 

N. Kidonakis: 

We were interested in the relative size of the NLL. Of course the result 
depends somewhat on the procedure used and the value of the cutoff. How- 
ever, we calculated the cross section with different cutoffs and found that the 
NLL contribution is always significant, and that is our main conclusion. 



